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ABSTRACT 

In a radio search with the Green Bank Telescope of three unidentified low Galactic latitude Fermi - LAT 
sources, we have discovered the middle-aged pulsar 12030+3641, associated with 1FGL J2030.0+3641 
(2FGL J2030.0+3640), Following the detection of gamma-ray pulsations using a radio ephemeris, we have 
obtained a phase-coherent timing solution based on gamma-ray and radio pulse arrival times that spans the 
entire Fermi mission. With a rotation period of 0.2s, spin-down luminosity of 3 x 10 34 ergs -1 , and character- 
istic age of 0.5 Myr, PSR J2030+3641 is a middle-aged neutron star with spin parameters similar to those of 
the exceedingly gamma-ray-bright and radio-undetected Geminga. Its gamma-ray flux is 1% that of Geminga, 
primarily because of its much larger distance, as suggested by the large integrated column density of free elec- 
trons, DM = 246 pc cm -3 . We fit the gamma-ray light curve, along with limited radio polarimetric constraints, 
to four geometrical models of magnetospheric emission, and while none of the fits have high significance 
some are encouraging and suggest that further refinements of these models may be worthwhile. We argue 
that not many more non-miiiisecond radio pulsars may be detected along the Galactic plane that are responsi- 
ble for LAT sources, but that modified methods to search for gamma-ray pulsations should be productive — 

PSR J2030+3641 would have been found blindly in gamma rays if only > 0.8 GeV photons had been consid- 
ered, owing to its relatively flat spectrum and location in a region of high soft background. 

Subject headings: gamma rays: observations — pulsars: individual (PSR J2030+3641) 


1. INTRODUCTION 

Three years after the launch of the Fermi Gamma-ray Space 
Telescope , pulsations have been detected with its Large Area 
Telescope (LAT) from at least 100 rotation-powered neutron 
stars (e.g., Abdo et al. 2010d), increasing by more than an 
order of magnitude the number of pulsars known to emit 
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above 0.1 GeV l7 . Nearly 20 of these are millisecond pulsars 
(MSPs) known independent from Fermi , for which rotational 
ephemerides obtained from radio observations were used to 
fold the sparse gamma-ray photons (e.g., Abdo et al. 2009a, 

2010a). Another 34 have been discovered in direct pulsa- 
tion searches of the gamma-ray data (e.g., Abdo et al. 2009b; 

Pletsch et al. 201 1 ; Saz Parkinson et al. 2010). 

Among the known non-MSP gamma-ray pulsar sample, 
about 75% have large values of spin-down luminosity ( E > 
10 35 ergs -1 ); the remaining 17 have lower E and are older, 
with characteristic age r c = P/(2P) ~ 0. 1-1 Myr. Five of these 
are also radio emitters (two are exceedingly faint and were 
detected only after the gamma-ray discoveries; Camilo et al. 

2009b; Pletsch et al. 2011). Almost all of these middle-aged 
pulsars are also known or expected to be nearby, at d < 1 kpc, 
although obviously many more similar pulsars exist at larger 
distances. In principle several of them could be detected in 
gamma rays despite a relatively low E , because the efficiency 
for conversion of rotational kinetic energy into gamma rays, 
t}=L 1 /E , apparently increases with decreasing E and may 
approach 100% for E « 10 33-34 erg s _1 (see Abdo et al. 2010d; 

Arons 1996; Muslimov & Harding 2003). 

In order to improve our understanding of pulsar emis- 
sion mechanisms and evolution, it is important to iden- 
tify more gamma-ray-emitting neutron stars, particularly in 
under-represented groups such as middle-aged pulsars. The 
Fermi LAT First Source Catalog (1FGL), based on 1 1 months 
of survey data and containing more than 600 unidentified 
sources (Abdo et al. 2010b), provides a path toward such 
discoveries (as does now 2FGL; Abdo et al. 2011a). Many 

17 See https://confluence.sIac.stanford,edu/display/GLAMCOG/Public+List+of+LAT- 
Detected+Ga mm a-Ray+Pulsars for an up-to-date list. 
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unidentified LAT sources are being surveyed within the con- 
text of & “pulsar search consortium” that aims to efficiently or- 
ganize this collaborative radio and gamma-ray work (see, e.g., 
Ray & Parkinson 201 1). Spectacularly, many radio MSPs 
have recently been discovered in such searches at high Galac- 
tic latitudes (|6| > 5°; e.g., Cognard et al. 2011; Keith et al. 
2011; Ransom et al. 2011). Here we report on the first pulsar 
to be discovered at low Galactic latitude that is responsible 
for a formerly unidentified 1FGL source. PSR J2030+3641 
in 1FGL J2030.0+3641 (also 2FGL J2030. 0+3640), discov- 
ered with the NRAO Green Bank Telescope (GBT), is also the 
first non-MSP gamma-ray pulsar identified in radio searches 
of 1FGL sources. 

2. OBSERVATIONS AND RESULTS 

2.1. Radio Searches 

Consistent with the properties of known gamma-ray pul- 
sars, our first 1FGL search targets are non- variable and have 
spectra consistent with exponentially-cutoff power laws. In 
this initial small radio survey at the GBT we aimed to search 
along tne Galactic plane, indicating a relatively high search 
frequency (> 1 GHz) in order to minimize deleterious propa- 
gation effects in the interstellar medium and the background 
temperature due to Galactic synchrotron emission. At the 
2 GHz frequency used, the GBT beam is 3' HWHM, and the 
1FGL sources to be searched should therefore, ideally, have 
rgs (95% confidence level error radii) below that, after ac- 
counting for the combined statistical and estimated system- 
atic positional uncertainties. These are strict criteria, and we 
searched a total of only three sources. For each we provide 
(a, <5, A^,r 9 5 ), which are respectively the R.A. and decl. of 
the observed position, its offset from the 1FGL position (we 
observed in late 2009, based on an interim version of the cata- 
log), and rgs for the gamma-ray source: 1FGL J0224.0+6201c 
(36P00,62?04, 1',4'X where the “c” indicates that the mea- 
sured properties of this source (including position) may not 
be reliable; IFGLJ 1746.7-3233 (266.°70,-32P61,3 , ,3 / ); and 
1 FGL J2030.0+364 1 (307.°52,36. o 68, l',3'). 

We observed each source for 1 hr using the GUPPI spec- 
trometer 18 at a central frequency of 2 GHz. Each of 512 
polarization-summed frequency channels, spanning a band- 
width of 800 MHz, were sampled every 0. 164 ms before writ- 
ing to disk. The flux density limit of these searches, converted 
to a standard search frequency of 1 .4 GHz, was approximately 
0.03 mJy, for an assumed pulsar duty cycle of 10% and spin 
period larger than a few tens of milliseconds. We analyzed 
the data with standard pulsar search techniques implemented 
in PRESTO (Ransom 2001), selecting trial dispersion mea- 
sures so as to maintain ideal sensitivity up to twice the maxi- 
mum Galactic DM predicted in each direction by the Cordes 
& Lazio (2002) electron distribution model. The data sets 
retained significant sensitivity to (possibly binary) MSPs up 
to DM « 100 pc cm -3 . We therefore did the analysis in two 
passes, including one where we allowed for a modest amount 
of unknown constant acceleration (parameterized in PRESTO 
by zmax = 50). In the 1FGL J2030.0+3641 data, collected on 
2009 November 27, we discovered on Christmas day an obvi- 
ous pulsar with period P = 0.200 s and DM = 247 pc cm 3 . 

2.2. Radio Timing and Polarimetry of PSR J2030+364 / 

1 8 https ://wikio.nrao.edu/bin/view/CIC ADA/GUPPiUsersGuide 



FIG. 1. — Polarimetric observation of PSR J2030+3641 at 1 .5 GHz with the 
GBT. In the lower panel the total-intensity profile is represented by the black 
trace, red corresponds to linear polarization, and blue to circular. The profiles, 
based on a 4.5 hr integration, are displayed with a resolution of 256 bins and 
arbitrary phase. In the upper panel the position angle of linear polarization 
is displayed for bins in which the linear signal-to-noise ratio exceeds 2.2, 
and the PAs have been rotated to the pulsar frame using the measured RM = 
+514 rad m -2 . The full pulse phase corresponds to P = 0.200 s. 

We began regular timing observations of the new pulsar at 
GBT on 2010 January 8. Observing parameters were identical 
to those used in the search observation, although with much 
shorter integration times, typically 5 minutes. In this man- 
ner we measured pulse times of arrival (TOAs) on 28 sepa- 
rate days through 201 1 April 1 1 . Using these with TEMPO 19 
we obtain a phase-connected rotational ephemeris with 0.3 ms 
rms residual. After 5 months of radio timing, the solution was 
sufficient to yield an initial detection of gamma-ray pulsations 
(see Section 2.3). The DM measurement was improved with 
the aid of one additional observation centered at 1 .5 GHz, and 
implies an uncertainty in the alignment of radio and gamma- 
ray pulses of only 0.7 ms. 

In order to determine the polarization characteristics of 
PSR J2030+3641, we have used GUPPI to make three flux- 
calibrated fiill-Stokes observations of the pulsar, one each at 
central frequencies of 0.82 GHz (lasting for 2.3 hr), 1.5 GHz 
(4.5 hr), and 2 GHz (l.Ohr). We analyzed the data with 
PSRCHTVE (Hotan et al. 2004), and show the resulting 
1 .5 GHz profiles in Figure 1 . The profiles at the other two 
frequencies are comparable, with a possible slight hint of 
interstellar scattering visible at 0.8 GHz. The Faraday rota- 
tion measure is RM = +5 14 rad nr 2 , which implies an average 
electron-weighted Galactic magnetic field strength along the 
line of sight of 2.6 /iG. The radio spectral index, based on 
the three flux density measurements, is o = -1 .7 ± 0. 1 (where 
S u oc v a \ see Table 1). 

2.3. Gamma-ray Study of PSR J2030+3641 

We used “Pass 6 diffuse”-class Fermi LAT events (having 
the highest probability of being gamma-ray photons), exclud- 
ing those with zenith angles > 100° to minimize gamma-ray 
emission from the Earth’s limb. The events were analyzed 
using the LAT Science Tools 20 , and photon phases were cal- 
culated using TEMP02 (Hobbs et al. 2006) with the f ermi 
plugin (Ray et al. 201 1). 

19 http://www.atnf.csiro.au/research/pulsar/icmpo 

20 http://fermi.gsfc.nasa.gov/ssc/data/analysis/scitools/overview.html 
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FIG. 2. — Phase-aligned radio and gamma-ray profiles of PSR J203O+3641. 
Bottom panel: discovery observation at the GBT. Middle panel: Fermi LAT 
0.5-0.8 GeV profile (significant LAT pulsations are not detected below this 
energy band). Top panel: > 0.8 GeV and > 3 GeV (lower) LAT profiles. 
Each pulse profile is repeated in phase, and displayed with 30 bins per period 
in gamma rays and 128 bins in radio. The dashed line gives the level of 
unpulsec emission in the aperture as estimated from the light curve fits. 

We first folded photons detected from 1FGL J203 0.0+3 641 
based on the initial radio timing solution for PSR J2030+3641 
(Section 2.2). Selecting photons using guesses for the mini- 
mum energy and radius of interest for extraction, we made 
a clear detection of pulsations. Optimal cuts in this re- 
gion of very high background were subsequently found with 
E > 0.8 GeV and within a 0.°5 radius. We then obtained 20 
gamma-ray TOAs from a total of 1402 photons detected be- 
tween 2008 August 4 and 2011 April 20 (49 days of integra- 
tion per TOA), and used TEMP02 to determine a timing solu- 
tion with an rms of 4.0 ms. Despite the much higher precision 
of the individual radio TOAs, the best overall ephemeris is 
obtained by combining the gamma-ray and radio TOAs, ow- 
ing to the longer span of the former. The parameters of this 
timing solution, in which a small amount of “timing noise” is 
parameterized by v (where v = 1 /P), are presented in Table L 

Figure 2 shows the gamma-ray pulse profiles in three en- 
ergy bands, along with the radio discovery profile. In order to 
optimize each pulsed signal-to-noise ratio, we selected pho- 
tons within 0.°8, 0?5, and 0.°3 of the timing position, respec- 
tively for the 0.5-0.8GeV, > 0.8 GeV, and > 3 GeV bands. 
For the middle energy band, the bin-independent H-test (de 
Jager & BUsching 2010), with a value of 270, indicates enor- 
mous significance. For each of the other two bands, the H- 
test values of approximately 40 correspond to pulsation sig- 
nificances over 5cr. The light curve at > 0.8 GeV consists of 
two (conceivably three) overlapping components as observed, 
e.g., for PSRs J1709-4429, J1057-5226, and the CTA 1 pul- 
sar (Abdo et al. 201 Od), which have peak separation A « 0.25. 
While the statistics are poor, the > 3 GeV profile is consistent 


with such a description. At 0.5-0. 8 GeV, only one peak is 
visible, which corresponds in phase with the first peak visi- 
ble at higher energies. No pulsed signal was detected below 
0.5 GeV, despite searches using different apertures. We fit 
the unbinned gamma-ray data above 0.8 GeV with two two- 
sided Gaussians, to account for the steep leading and trail- 
ing edges of the pulsed emission. The first peak is offset 
in phase from the main radio peak by <5 = 0.26 ± 0.02 with 
FWHM = 0.1 1 ± 0.04, while the second peak is offset from 
the first by A = 0.30±g;g|, with FWHM = 0.18 ± 0.05. We 
have included systematic error from the modest uncertainty 
in DM and our choice of model. 

2.3.1. Spectrum 

To characterize the spectrum of PSR J2030+3641, we an- 
alyzed Pass 6 diffuse events collected between 2008 August 
4 and 2011 April 20 selected to have reconstructed energies 
between 178MeV and 100 GeV and reconstructed positions 
within 10° of the timing position of PSR J2030+3641. We 
excluded the 1 00- 178 Me V energy bin because the rapidly 
increasing effective area in this bin exacerbates the effects of 
energy dispersion (Atwood et al. 2009). We removed events 
with a reconstructed zenith angle > 100° and also filtered data 
in which the aperture approached the limb of the Earth; these 
cuts remove most of the contribution of the bright gamma-ray 
emission produced by cosmic rays interacting in the Earth’s 
atmosphere. To further enhance the signal-to-noise ratio of 
PSR J2030+3641, we selected events with rotation phase 
0.22 < <b < 0.64 (see Figure 2). We verified that there was 
negligible emission outside of this phase range, so the best-fit 
values may be taken to represent the phase-averaged magne- 
tospheric emission. 

PSR J2030+3641 lies in the source-crowded Cygnus re- 
gion, and extracting its spectrum requires careful model- 
ing of both neighboring point sources (including the three 
bright gamma-ray pulsars J2021+3651, J2021+4026, and 
J2032+4127) and the diffuse background. We included all 
point sources in a preliminary version of the 2FGL catalog 
(Abdo et al. 2011a) within 15° of the pulsar, as the broad 
point-spread function (PSF) at low energy allows sources 
away from the aperture boundary to contribute to the observed 
counts. Identified pulsars were modeled with a power law 
with exponential cutoff, dN/dE oc (E / GeV) _r exp(-£ /E c ). 
In the fits reported below, we varied all parameters for sources 
within 7° of the aperture, while sources outside of the aperture 
were held fixed at the catalog values. We modeled the Galac- 
tic and isotropic diffuse backgrounds using, respectively, 
the gll_iem_v02_P6_Vl 1 -DIFFUSE diffuse map and the 
isotropic_iem_v02_P6_Vl 1 _DIFFUSE tabulation employed 
for the 1FGL spectral analysis; these have been rescaled from 
those constructed with the P6_V3_DIFFUSE instrument re- 
sponse functions (IRFs). We allowed the normalizations of 
the diffuse models to vary in the fit. 

We performed spectral fits with the pointlike application 
(Kerr 2011). This software bins the events in both energy and 
position, using HEALPix 21 to provide position bins whose 
energy-dependent extent is always small compared to the PSF, 
proviing a compact representation of the data with minimal 
information loss. We verified the results with the Science- 
Tool gtlike. The best-fit values are reported in Table 1, and 
include a power-law photon index T = 1 . 1 and cut-off energy 

21 http://healpix.jpl.nasa.gov/ 
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E c = 2.0 GeV. These values apply to the profile as a whole, 
although Figure 2 suggests that the second gamma-ray peak 
has a harder spectrum than the first, as is observed for many 
gamma-ray pulsars (e.g., Abdo et al. 2011b). 

In addition to the rescaled 1FGL diffuse model, we also per- 
formed the fit with an improved model internal to the LAT col- 
laboration and with “Pass 7” data 22 using the P7SOURCE_V6 
ERFs and the appropriate diffuse model (accompanying the 
2FGL catalog). These values were consistent with those re- 
ported in Table 1, and the scatter has been used to estimate 
systematic errors on the parameters. We have also included 
uncertainty in the instrument’s effective area in the system- 
atic errors through the use of “bracketing” ERFs (Abdo et al. 
2009c) 

2.4. X-ray Observations 

A short Swift (Gehrels et al. 2004) observation was ob- 
tained on 2010 May 1 1 to search for an X-ray counterpart to 
PSR J2030+3641. No emission is observed at the pulsar po- 
sition in this 3,8 ks X-ray Telescope (Burrows et al. 2005) PC 
mode observation, with a 3 u upper limit of 1.7 x 10 -3 cts s _1 
(0.5-8 keV). Assuming a power-law spectrum with photon 
index F = 1.5 for a possible pulsar wind nebula (e.g., Kar- 
galtsev & Pavlov 2008) absorbed by a column with Afe = 
7.4 x 10 21 cm -2 (assuming one free electron for every 10 neu- 
tral hydrogen atoms along the line of sight), this corresponds 
to an unabsorbed flux limit of fx < 1.3 x 10" 13 erg cm -2 s" 1 , or 
luminosity L x < 1.6 x lO 31 ^/! kpc) 2 ergs -1 . For a distance 
of « 2 <pc (see Table 1), this corresponds to L x /E < 2 x 10 -3 , 
which is a reasonably stringent limit, although it would not be 
a surprise if the actual X-ray flux of this pulsar were about an 
order of magnitude below the current limit (see, e.g., Marelli 
et al. 2011). 

3. DISCUSSION 

The gamma-ray pulsations from PSR J2030+3 641 unam- 
biguously identify it as the origin of the gamma rays from 
1FGL J2030.0+3641. We note also that a candidate TeV 
source possibly detected by the Milagro Gamma-Ray Ob- 
servatory (“C2” in Abdo et al. 2007) is positionally coinci- 
dent with PSR J2030+3641. Its faintness and location near a 
very bright source precludes further investigation based only 
on published results (see also Abdo et al. 2009d), but the 
relatively small E and relatively large age and distance of 
PSR J2030+3641 make an association unlikely. With a DM- 
derivec distance of 8kpc, PSR J2030+3641 is located toward 
the Cygnus region, where until recently very few pulsars were 
known and where the pulsar distance scale is highly uncer- 
tain. Recent discoveries suggest that the Cordes & Lazio 
(2002) model greatly overpredicts distances in this direc- 
tion: PSR J2032+4127, located 5° from PSR J2030+3641, 
with a DM-derived distance of 3.6 kpc, is thought to be 
locatec at half that distance (Camilo et al. 2009b); and 
PSR J2021+3651, less than 2° from PSR J2030+3641, with 
a model distance of 12 kpc, is more likely located at 2-4 kpc 
(Abdo et al. 2009e). Also, the RM for PSR J2030+3641 is 
large and positive, very similar to that of PSR J2021+3651, 
and this is suggestive of a location closer than the Perseus arm, 
at 6 kpc in this direction (see discussion in Han et al. 2006). 
Overall, we judge that a likely distance for PSR J2030+3641 
is in the range 2-4 kpc. 

22 See http://fermi.gsfc.nasa.gov/ssc/data/anaiysis/scitools/ 


3.1. Geometrical Constraints from Radio Profile 

Johnston & Weisberg (2006) noted that the profiles of pul- 
sars with high E share many common features. In particular, 
the profiles with two components generally show (i) two com- 
ponents with equal width, (ii) the trailing component stronger 
than the leading component, (iii) close to 100% linear polar- 
ization, (iv) significant circular polarization under the trail- 
ing component, and (v) flat position angle swings which ap- 
pear to steepen at the far trailing edge of the profile. The 
profile of PSR J2030+3641 conforms in some ways to these 
generalizations. The profile is symmetrical with the trailing 
edge brighter. The PA swing is remarkably flat with a hint of 
steepening at the trailing edge. However, although the leading 
component is highly polarized, the trailing edge is not, unlike 
the majority of the high E pulsars. Weltevrede & Johnston 
(2008) pointed out that the transition between pulsars with 
low linear polarization and those which are highly polarized 
occurs at E ~ 10 34 5 erg s -1 . PSR J2030+3641 could be an ex- 
ample of a pulsar in transition from high to low polarization. 

The interpretation of these profiles by Johnston & Weisberg 
(2006) was that the magnetic pole crossing occurs at the (sym- 
metry) center of the profile. The PA swing is significantly 
shifted to later phase because of aberration (Blaskiewicz et al. 
1991), leading to emission heights close to 1000 km. This 
may also be the case in PSR J2030+3641. If we assume that 
the inflexion point of the PA swing is the rise in the PA curve 
some 20° after the profile center in Figure 1 , this implies an 
emission height of ~ 800 km (close to 10% of the light cylin- 
der), similar to those seen in other high E pulsars. We also 
note that a height of ~ 800 km would in turn imply an overall 
pulse width of some 50°, similar to the observed value. This 
further suggests that the inclination between the magnetic and 
rotation axes cannot be too far from orthogonal. 

3.2, Light Curve Fitting and Model Implications 

We can use the gamma-ray light curve and polarimetric 
measurements to test the predictions of magnetosphere beam- 
ing models. Such tests are most constraining when the mag- 
netic inclination angle a and the viewing angle C are well 
determined. The polarization data provide some geometri- 
cal constraints. Since the lower frequency data may be af- 
fected by interstellar scattering, our primary constraints are 
based on the 1.5 GHz data. We use PA values from the na- 
tive 512 observational phase bins with a pa <10° and phases 
-0. 1 < <j) < 0.05 (relative to the intensity maximum; see Fig- 
ure 1), and fit to a standard rotating vector model (RVM) po- 
larization sweep. The maximum sweep rate is d^/d^max 
= 1 .4 and one obtains good fits in a band of the a-£ plane, 
with a best fit x 2 = 37.5 for 55 degrees of freedom (suggesting 
the PA errors are somewhat overestimated). Figure 3 shows 
the RVM fit confidence regions in this plane, after correcting 
for the RVM sign convention problem (Everett & Weisberg 
2001). The phase of maximum sweep is not well constrained, 
but the best fits place it well after the radio peak at <j>ta 0.06. 
The fit implies a < 140°, with best fit values < 90°. If we si- 
multaneously fit the 1.5 GHz points and the -0.1 < <j> < 0.02 
0.8 GHz data, larger a > 90° are preferred, but this preference 
is controlled largely by the trailing edge of the sweep where 
scattering distortions may be present. 

Following Romani & Watters (2010), we next compare the 
gamma-ray light curve to three geometrical models of outer 
magnetosphere emission. The first, “Outer Gap” (OG; Cheng 
et al. 1986), model posits emission from vacuum dipole field 
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Fig. 3. — Surfaces of goodness of fit (* 3 ) to the E > 0.8 GeV light curve for the OG (left) and SL (right) outer magnetosphere beaming models. Dark colors 
indicate lower (better) values of the fit statistic. The Ax 2 = +l,+2,+3 contours for the RVM fit to the 1.5 GHz polarization data are shown in green. 



FlG. 4. — Model light curves (red) at the best-fit a and < consistent with 
the RVM PA constraints compared to the E : 0.8 GeV LAT light curve (black 
histograms). The offset <p indicates the phase of the closest approach of the 
magnetic dipole axis at the star surface. The SL curves have been offset by 
+80 counts/bin. 

lines radiating from the null charge surface to the light cylin- 
der at radius /?lc- In this model, the heuristic efficiency scal- 
ing 7} ft# (10 33 ergs _1 /£) 1 ' 2 is assumed. For PSR J2030+3641, 
the value rj = 0. 18 was used for the gap width w in the OG 
calculations. The second model, “Two Pole Caustic” (TPC; 


Dyks & Rudak 2003), has emission in the vacuum open zone 
from the neutron star surface out to 0.75 /?lc from the rota- 
tion axis or an altitude of Rue , whichever is lesser. More re- 
cently, numerical models of force-free magnetospheres have 
been computed by Bai & Spitkovsky (2010) and they have ar- 
gued that the pulsar gamma-ray beam can be created in the 
“Separatrix Layer” where the open zone current sheets over- 
lap. The emission in this SL model is assumed to extend to 
1.5 Rlc- We assess the match to the data by plotting a good- 
ness of fit X 3 surface in the a-£ plane (see Romani & Wat- 
ters 2010). Using the E > 0.8 GeV light curve, the best-fit 
OG models (with best X 3 = 3.9 at a = 133°, £ = 1 10°) lie in- 
side the 3 a 1.5 GHz fit contours and at the maximum of the 
combined 1 .5 GHz and 0.8 GHz fit. A small second region of 
adequate but poorer fits (x 3 = 7) is also present within the 1 a 
RVM contour at smaller a = 63° (£ = 60°). The SL model 
has a similar small region of good fits within the 1 a 1 .5 GHz 
RVM contour (X 3 = 5 at a - 50°, £ = 48°). The TPC model 
predicts a large unpulsed component resulting in much poorer 
fits (not shown), with a best value X 3 = 19 in the RVM-allowed 
region; this is at a = 26°, £ = 25°, giving a light curve with 
little pulsed flux. 

Figure 4 shows the best-fit light curves in the RVM-allowed 
region for the OG and SL models, compared with the data. 
The present models do not include the full radiation physics, 
employing simple uniform emissivity throughout the active 
zones. However, it appears that both high altitude models can 
produce reasonable matches to the LAT data. With improved 
geometrical constraints it may be possible to distinguish be- 
tween the OG and SL models. For the former, the fit phase of 
the magnetic axis is at 0 ft# 0.09, after the pulse and somewhat 
later than the polarization sweep. The SL model prefers a 
magnetic axis much closer to the maximum of the radio pulse. 
In both cases, the lag of the sweep maximum from the pulse 
peak implies significant altitude ~ 0. 1-0. 15Rlc for the radio 
emission. Such altitudes are also required to accommodate 
the relatively wide radio pulse in the open zone. Modeling 
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a 

FIG. 5. — Reduced x 2 values (with 25 degrees of freedom) in a-£ space 
from slot gap gamma-ray light curve fits, marginalized over w , r, and A <t>, and 
black curves from the RVM fit to radio, polarization measurements, showing 
where these constraints intersect the geometrical light curve parameters (the 
overall best-fit light curve lies outside these regions; see Section 3.2). The 
plotted values of x 2 s ^ those included in the Markov chains via the algo- 
rithm given in Verde et al. (2003). 


such finite altitude effects may provide some additional re- 
strictions on the geometry. We also note that the OG and SL 
models predict somewhat different viewing angles (; if an X- 
ray pulsar wind torus could be measured, the resulting ( might 
be used to select the preferred model. 

The classic TPC model fails in these fits at least in part be- 
cause it does not have enough emission at high altitudes. In 
order to improve upon this, we also fit the > 0.8 GeV profile 
of PSR J2030+3641 with light curves predicted from a ge- 
ometrical high-altitude “Slot Gap” model (SG; Muslimov & 
Harding 2004), an evolution of the low-altitude SG of Mus- 
limov & Harding (2003). In the SG geometry, the emission 
originates in a gap that extends from the neutron star surface 
at the polar cap to high altitudes near the light cylinder. We fit 
for four model parameters; a, gap width w (a fraction of the 
open volume radius of the polar cap, r ovc , in the open volume 
coordinate system described in Dyks & Harding 2004), and 
maximum radius of emission r (in units of Rlc\ assuming 
uniform emission rate in the corotating frame. Our resolu- 
tion is 1° in ol and (, 0.01 r ovc in 0 < w < 0.3, and 0.1 R L q 
in 0.71?lc < r < 2.0 R^c (we assumed a maximum cylindrical 
emission radius of 0.95 The simulated light curves are 
fit to the data via a Markov Chain Monte Carlo (MCMC) rou- 
tine that explores the parameter space and records regions of 
maximum likelihood, as in Verde et al. (2003). 

As noted before, we have some independent constraints on 
the emission geometry of PSR J2030+3641 from radio polar- 
ization (cf. Figure 1). To make our model fits be consistent 
with the RVM fit, we allowed the model light curve to shift 
with respect to the location of the magnetic pole by no more 
than = 0.1 in phase. Figure 5 shows the a- £ space with 
the polarization contours already shown in Figure 3 and color 
“contours” of reduced x 2 from the MCMC likelihood fits. The 
polarization contours pass through parameters yielding light 
curve fits with 3 < x 2 /25 < 5, the best fit falling within the 
contours occurring at ( a , £, w, r, A 0; x 2 /25) = (134°, 1 14.5°, 
0.0 rovc* 1-4Rlc, 0.08; 3.4). The second best fit that falls 
within the polarization contours is at (54°, 52.5°, 0.01 r ovg , 
1.3/?lc> 0.105; 4.2). We include this fit because while it has 
a larger x 2 > it falls within a more constrained region of the 


polarization contours. The absolute best fit, with or without 
restricting A0, is at (170°, 110.5°, 0.0r ovc , > l.l/? L cs 0.017; 
1.4). This is inconsistent with standard radio cone beam emis- 
sion because the values of a and £ are too far apart for lower 
altitude radio emission to be visible. All of these light curves 
are shown in Figure 6. 

Two of the fits include w = 0, meaning the emission in the 
model originates from a single open field line. Physically, this 
means the emission zone width is smaller than the resolution 
of our simulations, but not really zero. All the fits have small 
widths, requiring higher current of radiating particles in the 
gap (either more primaries or high pair multiplicity) in order 
to generate the observed luminosity. It is also likely that the 
emissivity and width of the emission zone change azimuthally 
and with altitude, unlike the constant width of our simula- 
tions (see Harding & Muslimov 2011). None of the fits are 
able to reproduce well the small middle peak of the observed 
light curve; it is not yet clear whether this is a separate feature 
or part of the second peak. Also, the background is slightly 
over-predicted by the SG model, especially for the parameters 
consistent with the polarization curves. 

None of these (OG, SL, SG) gamma-ray light curve fits are 
particularly consistent with our inferences from the radio pro- 
file and polarimetric constraints. For instance, the latter sug- 
gest a nearly orthogonal rotator (Section 3.1), for which there 
are no good gamma-ray light curve fits. However, the po- 
larimetric data provide limited constraints, owing to the small 
range of pulse longitude for which they are obtained, and our 
inferences rely on some assumptions. As for the gamma-ray 
fits, none are perfect, although the OG and SG fits appear 
slightly better than the SL one. Considering their inherent 
limitations, however, these fits are encouraging, and further 
refinements may improve our understanding of the magneto- 
sphere of middle-aged pulsars like PSR J2030+3641. 

For each beaming model, one can extrapolate the aver- 
age pulsar flux on the Earth line-of-sight to the isotropic all- 
sky equivalent. For the best-fit angles consistent with polar- 
ization constraints, the corrections for PSR J2030+3641 are 
jfo = 0.66 (OG), 0.64 (SL), and 0.78 (SG; neglecting polar- 
ization constraints, /q,sg = 0.45). Using the observed flux 
and the 8kpc DM-estimated distance, we find that r) 0 bs = 
AitdtfnF^E = 10dg/h. This indicates that any distance 
greater that ~ 2.5 kpc is problematic, unless the emission is 
highly beamed. For the OG/SL/SG beaming fy ~ 0.7, we re- 
quire d < 3 kpc. The preferred distance for the heuristic OG 
efficiency law (77 = 0. 1 8) is ~ 1.3 kpc. Such distance estimates 
of « 2 kpc are consistent with the discussion at the beginning 
of Section 3, and strongly suggest that PSR J2030+3641 is lo- 
cated at roughly lOx the distance to Geminga ( d = 0.16 kpc). 
Improved distance estimates for middle-aged, apparently effi- 
cient pulsars have the potential to constrain or falsify gamma- 
ray emission models. 

3.3. Where Have All the ( Ordinary Radio) Pulsars Gone? 

In a radio search with the GBT of three unidentified 
low Galactic latitude LAT sources, we have discovered 
PSR J2030+3641, a relatively distant member of the under- 
represented group of gamma-ray pulsars that are middle- 
aged and have small E. A gamma-ray-only pulsar was sub- 
sequently discovered in one of the two remaining sources, 
1FGL J 1746.7-3233 (Pletsch et al. 2011). The third source 
could yet be a gamma-ray pulsar, but is unlikely to be a 
substantial radio emitter beamed towards the Earth (assum- 
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Frc. 6. — Slot gap model light curves compared to the data, (a) The best-fit light cun e (blue) with parameters lying within the polarization contours. The LAT 
light curve is shown in black. The vertical dot-dashed line shows the model location of the magnetic pole, shifted by A<p from zero. The horizontal dashed line 
represents the background count level, (b) Same as (a), for the second-best fit parameters (red), (c) Same as (a) and (b) for the absolute best-fit light curve (green) 
found by the MCMC routine, with no restriction on A <j>. Its parameters do not fall within the radio polarization contours. 


ing we covered its location adequately — in the case of 
1FGL J1746.7-3233, the pulsar is 3f8 offset from our GBT 
search position; see Section 2.1): for a fiducial distance of 
4 kpc, our sensitivity corresponds to a luminosity limit L\ A — 
S lA d 2 = 0.5 mJy kpc 2 , which is near the low end of the radio 
luminosity distribution for detected radio pulsars (see Abdo 
et al. 2010c; Camilo ct al. 2009a, b, which include the two ex- 
ceptions with much lower luminosities). It would seem that 
we were quite successful, with one related radio discovery 
among three LAT sources searched. However, as we argue 
next, we don’t expect such continued success. 

Many more Galactic plane unidentified LAT sources have 
been searched for radio pulsars than the three mentioned here. 
This work is being done at the GBT, Parkes, Nang ay, GMRT, 
and Effelsberg telescopes by members of the Fermi pulsar 
search consortium, and includes searches of all sources that 
have pulsar characteristics at all Galactic latitudes. So far this 
has resulted in the discovery of more than 30 MSPs, all at high 
latitudes (e.g., Cognard et al. 2011; Keith et al. 201 1 ; Ransom 
et al. 201 1) — and PSR J2030+3641 . While a detailed popu- 
lation analysis will have to wait for the publication of all the 
completed searches, the lone example of PSR J2030+3641 al- 
ready stands out. Why have more radio pulsars not been de- 
tected among the unidentified Galactic plane LAT sources? 

A point related to this question concerns the very deep ra- 
dio searches done of 34 pulsars discovered in blind pulsa- 
tion searches of gamma rays (e.g., Abdo et al. 2009b; Pletsch 
et al. 2011; Saz Parkinson et al. 2010): only four of them 
were found to also be radio emitters (see Pletsch et al. 2011; 
Ray et al. 2011). Three are exceptionally faint and their de- 
tection was possible only following the gamma-ray discov- 
eries (Abdo et al. 2010c; Camilo et al. 2009b; Pletsch et al. 
201 1). The remaining one, PSR J2032+4127 (Camilo et al. 
2009b), is like PSR J2030+3641 in that they are both located 


in the Cygnus region and both are relatively bright, being de- 
tectable in 1 minute with the GBT at 2 GHz. They were not 
known prior to the launch of Fermi simply because no sensi- 
tive surveys had (or have yet fully) been done of the northern 
(equatorial) reaches of the Galactic plane, using the high radio 
frequencies desired to unveil the Galactic disk. Conversely, 
the enormously successful 1 .4 GHz Parkes multibeam survey 
(e.g., Manchester et al. 2001) did search the entire Galactic 
plane at 260° < / < 50° with sensitivity sufficient to detect or- 
dinary (non-MSP) radio pulsars with L\ A > 0.5 mJy kpc 2 out 
to d > 2 kpc. A very large fraction of all potentially detectable 
ordinary radio pulsars in that area of the sky that could also be 
plausible gamma-ray sources were therefore discovered be- 
fore Fermi launch, especially when adding those discovered 
in deep pointed observations of known pulsar wind nebulae 
(e.g., Camilo et al. 2002). In time, several of these pulsars 
were detected by LAT (e.g., Abdo et al. 2010d). Subsequent 
to the Parkes multibeam survey an even deeper Parkes survey 
extended Galactic plane coverage out to / = 60° (Camilo et 
al., in preparation), and the on-going PALFA Arecibo survey 
is doing so out to l « 75° (Cordes et al. 2006). 

In our view, the above suggests that the answer to “why 
have more radio pulsars not been detected among the uniden- 
tified Galactic plane LAT sources?” is that the vast majority 
of ordinary radio pulsars accessible to the current generation 
of telescopes and located within a few degrees of the plane at 
/ < 75° that can yield an appreciable gamma-ray flux at the 
Earth were discovered a long time ago (it is possible that in 
rare but important cases very high Ly distant pulsars will have 
their radio pulses scattered beyond practical detectability). 

Many of the LAT sources that remain unidentified along 
the Galactic plane are surely pulsars — but they may not be 
detectable as radio sources, and should be searched anew in 
gamma rays (some of these may be in binaries, and inacces- 
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sible tc current blind searches; see, e.g., Corbet et al. 2011). 
PSR J2 030+3641 is bright enough that it would have been 
discovered in blind searches of 18 months of LAT data, if 
the photon selection criteria had been adjusted to take into ac- 
count its flat spectrum (e.g., if only photons with E > 0.8 GeV 
had been searched for pulsations), as we confirmed after dis- 
covery. Spectral analysis of a region can improve the signal- 
to-noise ratio by selecting events based on the probability that 
they come from the source of interest (see Belfiore 2011), and 
applying this technique to the Cygnus region also resulted in 
the unbiased detection of pulsations from PSR J2030+3641. 
For fainter point sources superimposed on the large and un- 
certain diffuse Galactic background, spectral analysis and lo- 
calization are harder (1FGL J2030.0+3641 was only 0f4 from 
the actual pulsar position), but these and other improvements 
in search techniques are already resulting in the discovery 
of many more gamma- ray pulsars (see Pletsch et al. 2011; 
Saz Parkinson 2011). As a consequence, the ratio of known 
gamma-ray-only to gamma-ray-and-radio ordinary pulsars, 
which is currently slightly under 1.0, should increase substan- 
tially. 

When will we ever learn all that we can (in the Fermi 
era) about the geometry and emission properties of gamma- 
ray and radio pulsar accelerators? Perhaps when full details 
emerge from the continuing collaborative radio and gamma- 
ray observational and modeling work. The future is bright (at 
some wavelength), but requires publication of all searches, 
including radio non-detections, and consideration of the non- 
detection by LAT of radio pulsars with very large E flux at the 
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TABLE 1 

Measured and Derived Parameters for PSR J2030+3641 


Parameter 


Value 


Right ascension, R. A. (J2000.0) 

Declination, decl. (J2000.0) 

Galactic longitude, l (deg) 

Galactic latitude, b (deg) 

Rotation frequency, v (j _1 ) 

Frequency derivative, v (s“ 2 ) 

Frequency second derivative, 0 (j“ 3 ) 

Epoch of frequency (MJD) 

Dispersion measure, DM (pc cm -3 ) 

Data span (MJD) 

RMS timing residual (ms) 

Spin-down luminosity, E (ergs -1 ) 

Characteristic age, t c (yr) 

Surface dipole magnetic field strength (Gauss) — 

Flux density at 2 GHz, $2 (mJy) 

Flux density at 1.5 GHz, S 1.5 (mJy) 

Flux density at 0.8 GHz, S 0.8 (mJy) 

Rotation measure, RM (rad nr 2 ) 

Radio-gamma-ray profile offset, 6 (P) 

Gamma-ray profile peak-to-peak separation, A {P) 

Gamma-ray (> 0. 1 GeV) photon index, T 

Gamma-ray cut-off energy, E c (GeV) 

Photon flux (> 0.1 GeV) (lO^cm^s 1 ) . . ! 

Energy flux (: 0. 1 GeV), F 1 (10' H erg cm -2 s“ l ) . 

DM-derived distance (kpc) 

Most likely distance®, d (kpc) 


20 h 30 m 00?261(4) 

36°41'27"15(7) 

76.12 

-1.44 


4.99678736138(4) 

-1.62294(5) X 10“ 13 

2 .2(6) x 10" 24 

55400.0 

246.0(7) 

54700-55671 

0.4 


3.2 x 10 34 
4.9 x 10 s 

1.2 x 10 12 
0.09 
0.15 
0.40 


+514 ±1 
0.26 ±0.02 


0.30^81 

l.l±0.2^j 

2.0 ± 0.3^4 
3.6±0.5±0.8 


4.2 ±0.3 ±0.5 
8 

1,5-3 


NOTE. — Numbers in parentheses represent the nominal 1 <r TEMP02 
timing uncertainties on the last digits quoted. For gamma-ray parameters, 
the first uncertainty is statistical and the second accounts for systematics (see 
Section 2.3.1). 

a See discussion at the beginning of Section 3 and at the end of Section 3.2. 





